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Abstract. 

We show how models of neutrino masses and mixings can be differentiated on the basis of their predictions for 6*13 and 
lepton flavor violation in radiative charged lepton decays and fl — e conversion. We illustrate the lepton flavor violation results 
for five predictive 50(10) SUSY GUT models and point out the relative importance of their heavy right-handed neutrino mass 
spectra and 613 predictions. 
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1. INTRODUCTION 

Many models exist in the literature which attempt to explain the observed neutrino masses and mixings. The viable 
models agree with the presently known neutrino oscillation parameters falling within the 2(7 ranges flJJ] : 

sin 2 0i2 = 0.28 - 0.37, Am^ = (7.3 - 8.1) x 10~ 5 eV 2 , 

sin 2 13 = 0.38 - 0.63, Am^ = (2.0 - 2.8) x 10~ 3 eV 2 , (1) 
sin 2 0i3 < 0.033. 

The data suggests the approximate tri-bimaximal mixing texture of Harrison, Perkins, and Scott J2]: 

2/V6 1/V3 



Upmns=\ -1/V6 1/V3 1/V2 ), (2) 
-1/V6 I/a/3 1/V2 



with sin 2 23 = 0.5, sin 2 [2 = 0.33, and sin 2 13 = 0. 

The reason for the plethera of models still in agreement with experiment of course can be traced to the inaccuracy 
of the present data. In addition, there are a number of unknowns that must be still determined: the hierarchy and 
absolute mass scales of the light neutrinos; the Dirac or Majorana nature of the neutrinos; the CP-violating phases 
of the mixing matrix; any departure of the reactor neutrino angle, 0i 3 from 0°; any departure of the atmospheric 
neutrino mixing angle, 023, from 45°; whether the approximate tri-bimaximal mixing is a softly-broken or accidental 
symmmetry; and the magnitude of the charged lepton flavor violation? In this presentation we survey some of the 
models to determine what they predict for 0i 3 , the hierarchy, and lepton flavor violation. 



2. THEORETICAL FRAMEWORK AND MODELS 

The observation of neutrino oscillations implies that neutrinos have mass, with the mass squared differences given in 
Eq.(Q]). Information concerning the absolute neutrino mass scale has been determined by the combined WMAP, SDSS, 
and Lyman alpha data [3J which place an upper limit on the sum of the masses, < 0.17 eV. An extension of the 
SM is then required, and possible suggestions include one or more of the following: 

• the introduction of dim-5 effective non-renormalizable operators; 
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• the addition of right-handed neutrinos with their Yukawa couplings to the left-handed neutrinos; 

• the addition of direct mass terms with right-handed Majorana couplings; 

• the addition of a Higgs triplet with left-handed Majorana couplings. 

The general 6x6 neutrino mass matrix in the B(v a L, N^ L ) flavor basis of the six left-handed fields then has the 
following structure in terms of 3 x 3 submatrices: 



M L M T N 

M N Mr 



(3) 



where Mn is the Dirac neutrino mass matrix, Ml the left-handed and Mr the right-handed Majorana neutrino mass 
matrices. With Mi = and Mn « Mr, the type I seesaw formula is obtained |4(] 

M v = -M T N M R l M N , (4) 

for the light Majorana neutrinos, while if Ml ^ and Mn <<Mr, one obtains the type II seesaw formula Q5I], 

M v =M l -MIm r x M n . (5) 

The effective light Majorana mass matrix is complex symmetric and can be diagonalized by a unitary transformation, 
U vL , to give 

M^'" 8 = Uy L M v U V L = diag(mi , m 2l m 3 ), (6) 

with real, positive masses down the diagonal. On the other hand, the Dirac charged lepton mass matrix is diagonalized 
by a bi-unitary transformation according to 

M E ag = uIrM e Uel = diag(m e , m^m T ). (7) 
The neutrino mixing matrix, Vpmns [6], is then given by 

VpMNS = Ut L U v i — UpMNS^i 

(8) 

4> = diag(e'*i, e ! * 2 ,l), 

in terms of the approximately tri-bimaximal mixing matrix, Upmns an d the phase matrix, <t>, since an arbitrary phase 
rotation of U v l is not possible in the above. 

Models which have been introduced to explain the neutrino oscillation phenomena fall into two categories. There are 
those based on some lepton flavor symmetry such as fi — T interchange symmetry, the more restrictive S 3 or A4 flavor 
symmetry, and SO(3) or SU (3) flavor symmetries. Many attempts have been made to explore the location of texture 
zeros for the lepton mass matrices, with the hope that some flavor symmetry might be identified that way. A more 
ambitious class of models is based on SU(5) or 5(9(10) grand unification, where one attempts to explain the masses and 
mixings in the quark sector as well as the lepton sector. These models are said to have a "minimal" Higgs structure, if 
the Higgs bosons responsible for electroweak symmetry breaking transform as the 10, 126 dimensional representations 
of 5(9(10), and possibly the 120, 45 or 54 representations. They lead to symmetric or antisymmetric entries in the 
mass matrices. Other 5(9(10) models have Higgs bosons which transform as 10, 16, 16, and 45 representations and 
are referred to as "lopsided," since lopsided contributions to the down quark and charged lepton mass matrices can 
occur due to the SU (5) structure of the 16's. For recent reviews, cf. Ref. [7]. 



3. SURVEY OF MODEL PREDICTIONS FOR 613 

In a previous publication [8] we have made a survey of 63 models in the literature which give the large mixing 
angle (LMA) solution for the solar neutrino oscillations and firm and reasonably restrictive predictions for the reactor 
neutrino angle. The cutoff date for the selection of models was chosen to be May 2006. In this study we found most of 
the models predict a value in the range of 10~ 4 < sin 0i3 < 0.04 with a normal hierarchy preferred by 3: 1 . The results 
are displayed in Fig. 1 in the form of a sin 2 613 histogram, where the models are distinguished according to their type, 
and each is assigned the same area on the histogram though their imprecise predictions may cover several bins. 
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FIGURE 1. Histogram of sin 613 predictions for all 63 models. 



Since the planned Double CHOOZ and Daya Bay reactor experiments [9] will reach sin 2 20i3 ~ 0.01, roughly half 
of the models will be eliminated, either by the positive or negative signal of disappearance of v e 's from their beams. 
But even if a positive signal for the disappearance is seen, Fig. 1. indicates that of the order of 5 - 10 models will still 
survive, depending upon the accuracy of the measurement and the resilience of the models. 

Meanwhile the MEG experime nt [1 1 CHI at PSI is beginning to look for the \i — > <?y decay mode. With plans to lower 
the present branching ratio limit 111 ill of 1.2 x 10~ n down to the 10~ 13 range, this experiment may serve as an even 
more immediate selector of models. With this in mind, we turn to the subject of charged lepton flavor violation as a 
further distinguishing feature of the models proposed. 




FIGURE 2. Example of a Feynman diagram for fl 



ey with a neutrino mass insertion in the SM. 



4. LEPTON FLAVOR VIOLATION IN RADIATIVE DECAYS 



It is of interest to look first at charged lepton flavor violation in the SM with the addition of three massive right-handed 
neutrinos. In this case, individual L e , and L T lepton numbers are not conserved since Majorana mass terms are 
not forbidden. The lepton flavor violation then arises in 1-loop diagrams, where the neutrino insertion involves lepton 
flavor-changing Yukawa couplings. Examples of such diagrams are given in Fig. 2 for the process pi — > ey. The t/ft are 
elements of the PMNS mixing matrix, while the are Yukawa couplings for i,k= 1,2,3, and Nf is one of the heavy 
left-handed conjugate neutrinos. The branching ratio is given by 11211 



BR21 



r(n -> ey)/V{n -» v^eVe) 



3a 
31% 



3al28 



2 \2 



(9) 



sin 2 20 



12 ' 



10 



-54 



In this extended SM, the branching ratio is immeasurably small, due to the approximate GIM cancellation and the 
extremely small mass ratios of the left-handed neutrinos to the massive right-handed neutrinos. For such models, the 
MEG experiment would be expected to give a null result. 

In SUSY GUT models, slepton - neutralino and sneutrino - chargino loops contribute to the radiative decays as 
shown in Fig. 3. 
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FIGURE 3. Examples of Feynman diagrams for slepton - neutralino and sneutrino - chargino contributions to fl — > ey in SUSY 
models with slepton mass insertions. 

In the constrained MSSM (CMSSM) with universal soft masses and trilinear couplings lfl3ll . charged lepton flavor vio- 
lation arises from evolution of the Yukawa couplings and soft SUSY-breaking parameters. Thus with more comparable 
heavy masses in the loops and no GIM mechanism anticipated, the LFV branching ratios can be much larger [ 14]. 

In the leading log approximation the largest contribution comes from the left-handed slepton mass matrix yielding 
a branching ratio 



BR21 



G 2 F ml 



' 1 LL)J'\ 



Wj3, 



where 



( m LL)ji = ~^2 OT o( 3 



-Al/rnl)Yj k lo g m 



Y ki . 



(10) 



(ID 



with the Yukawa couplings specified in the lepton flavor basis and the right-handed Majorana matrix diagonal, so 
is just the kth heavy right-handed neutrino mass, while Mq is the GUT scale typically equal to 2 x 10 16 GeV and m s 
is some typical SUSY scalar mass. Petcov and collaborators [15] have shown that the full evolution effects as first 
calculated by Hisano, Moroi, Tobe, and Yamaguchi IU6I1 can be extremely well approximated by Eq. ( [Tol l, if one sets 



.0.5mlM 2 1/2 (m 2 - 



-0.6M? 



1/2J 



(12) 



We shall see that for the SUSY GUT models to be considered, the MEG experiment will be able to observe the 
predicted jj. — > branching ratios or place further restrictions on those models. 



5. GENERIC APPROACH TO LEPTON FLAVOR VIOLATION IN SUSY GUTS 



Many papers have studied predictions for the lepton flavor violating processes in SUSY GUT models by adopting a 
generic approach ifnll . Following the procedure of Casas and Ibarra lfl8ll . with the charged lepton and right-handed 
Majorana neutrino mass matrix diagonal, the seesaw formula can be inverted to yield the Yukawa neutrino coupling 
matrix 

Y v = —^DMiVM^RDyi^U^ (13) 
vsinp 

in terms of a complex orthogonal R matrix which allows for various unknown right-handed neutrino mixings. By using 
soft SUSY-breaking benchmarks, adopting various heavy right-handed neutrino masses and various R parametric 
angles and phases, one can "predict" the radiative LFV branching ratios. The results are typically presented in the 
form of Monte Carlo scatter plots. One finds the results strongly depend on tan j3 , M3 , 613, and the R parameters. The 
branching ratios are found to increase by powers of 10 as 613 varies from 0° — > 10°, while the present branching ratio 
bound for jU — > ey restricts M 3 to lie lower than 10 14 — 10 15 GeV. 

6. EXAMPLES OF PREDICTIVE SUSY GUT MODELS 

Instead of adopting a generic approach as described above, we are interested in determining whether one can differen- 
tiate between various SUSY GUT models on the basis of their LFV predictions, even if they have similar predictions 
for sin 2 613. For this purpose we have selected five 5(9(10) SUSY GUT models which are highly predictive, i.e., their 
model parameters are precisely specified by the authors. 

The models differ in the flavor symmetry chosen and the Higgs representations used to break the SO (10) symmetry 
at the GUT scale and the electroweak symmetry at the weak scale. We list the models, references, flavor symmetry 
and Higgs representations. 

• AB (Albright-Barr) H]: U(l) xZ 2 xZ 2 with 10, 16, 16, 45 

• CM (Chen-Mahanthappa) & SU(2) x (Z 2 ) 3 with 10, 126 

• CY (Cai-Yu) [21]: S A with 10,126 

• DR (Dermisek-Raby) JH]: D 3 with 10, 45 

• GK (Grimus-Kuhblok) [23]: Z 2 with 10, 120, 126 

Salient features of each model are listed in Table I. The CM model has a relatively large prediction for sin 2 613 which 
should be easily accessible to the future Double CHOOZ and Daya Bay reactor experiments. The AB, CY and DR 
model predictions for sin 2 613 ~ 0.0025 are barely within reach of those experiments. The GK model prediction would 
most likely require a neutrino factory to reach that level. Note that relatively low values of tan/3 are preferred for 
four of the models, while the DR model favors a high value. Relatively mild hierarchies for the massive right-handed 
neutrinos are predicted for the CM, DR and GK models, but the heaviest one occurs in increasing order for these three 
models, ranging from 7 x 10 12 to 2 x 10 15 GeV. The CY model predicts a degenerate right-handed spectrum with a 
lower value of 2.4 x 10 12 . The AB model, on the other hand, puts the heaviest one at 2.4 x 10 14 GeV and the lower 
two nearly degenerate at 4.5 x 10 8 GeV. Resonant leptogenesis [24] is possibly an interesting feature of that model. 
The five predictive models studied in this work thus cover a wide range of possibilities and suggest that charged lepton 
flavor violation can play an important role in further narrowing the list of viable candidates. 

7. RADIATIVE LEPTON FLAVOR VIOLATION PREDICTIONS 

We now turn to the predictions for the branching ratios for ji — > ey, % — ► jiy and T — ► ey in the five predictive models 
under consideration. We use the shorthand convention BR21, BR32, and BR31 for these respective branching ratios. 
Working in the CMSSM scenario with universal soft parameters mo, M\h and Aq for a given tanj3 and sgn(ji), we 
find a variety of plots are possible: 

• BRij vs. M1/2 f° r fixed Ao = and different choices of mo. 

• Ao/mo vs. Mi /2 scatterplot with a color scheme to indicate the branching ratio ranges. 



TABLE 1. Higgs representations, flavor symmetries, and other noteworthy features of the five 50(10) SUSY GUT 
models considered in this work. 



Models 


Higgs Content 


Flavor Symmetry 


M R (GeV) 


tan/3 


sin 2 0i3 


Interesting Features 



AB 


10, 16, 16,45 


U(\) xZ 2 xZ 2 


2.4 x 10 14 

4.5 x 10 8 
4.5 x 10 8 


5 


0.0020 
(2.6°) 


Large Mr hierarchy with lightest 
two nearly degenerate leads to 
resonant leptogenesis. 


CM 


10, 126 


5(7(2) x (Z 2 ) 3 


7.0 x 10 12 
4.5 x 10 9 

1.1 x 10 7 


10 


0.013 

(6.5°) 


Large Mr hierarchy with heaviest 
more than 3 orders of magnitude 
below GUT scale; large sin 2 613. 


CY 


10, 126 


Sa 


2.4 x 10 12 
2.4 x 10 12 
2.4 x 10 12 


10 


0.0029 
(3.1°) 


Degenerate Mr spectrum 4 orders 
of magnitude below GUT scale. 


DR 


10, 45 


D 3 


5.8 x 10 13 
9.3 x 10 11 
1.1 x 10 10 


50 


0.0024 
(2.8°) 


Mild Mr hierarchy almost 3 orders 
of magnitude below GUT scale. 


GK 


10, 120, 126 


z 2 


2.0 x 10 15 

4.1 x 10 14 
6.7 x 10 12 


10 


0.00031 
(1.0°) 


Mild Mr hierarchy just 1 order of 
magnitude below GUT scale; 
rather small sin 2 613. 



Ratio of the branching ratios, BR32/BR21 and BR31/BR21 on log-log plots, where for example, 

2 



logB/?32 = logfifl21 +log 



(i^v)32 



21 



(14) 



with L = log(M?/M|). Due to the factorization of the soft-breaking parameters and the GUT model parameters 
in the approximation of Eqs. (fTTT i and (fT2l . the above yields a straight line with unit slope and intercept given by 
the second term on the right. The length of the straight line depends on the soft parameter constraints applied. 

For the sake of brevity, we shall present only the third type of plots. For a more complete description we refer the 
reader to a longer paper submitted for publication i25ll . 

We have imposed the following soft parameter constraints lf26ll : 



For tan/3 =5,10: m : 50 -► 400 GeV 

M 1/2 : 200 -► 1000 GeV 
A Q : -4000 -► 4000 GeV 



Fortanj3=50: m Q : 500 -> 4000 GeV ( 5 

M 1/2 : 200 -> 1500 GeV 
A : -50 -> 50TeV 

In addition it is desirable to impose WMAP dark matter constraints [30 in the neutralino, stau or stop coannihilation 
regions l27ll . where the lightest neutralino is the LSR These more restrictive constraints are well described by the 
quadratic polynomial for the soft scalar mass in terms of the soft gaugino mass l28ll : 

m = c + ciM 1/2 + c 2 M 2 /2 , ^ 
where c; = c;(An, tan j3, sgn(ji)). 

where mo is bounded since M l / 2 is bounded. If M l / 2 is too small, the present experimental bound on the Higgs mass 
of rtih > 1 14 GeV may be violated or the neutralino relic density in the early universe will be too small, while if M^ 2 
is too large the neutralino relic density will be too large. 
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FIGURE 4. Branching ratio predictions for x — > fl + y vs. branching ratio predictions for jl — > e + y in the five models considered. 
The soft SUSY breaking constraints imposed apply for the thin line segments, while the more restrictive WMAP dark matter 
constraints apply for the thick line segments. The present experimental constraints are indicated by the dashed lines. 
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FIGURE 5. Branching ratio predictions for x — > e + y vs. branching ratio predictions for fl — > e + y in the five models considered. 
The same additional conventions apply as in Fig. 4, 



In Figs. 4 and 5 we have plotted BR32 and BR31 vs. BR21 on log-log graphs. The thin line segments for each 
model observe the soft parameters constraints imposed, while the heavier line segments observe the more restrictive 
WMAP dark matter constraints. The vertical dashed line reflects the present BR21 bound lUlll . while the horizontal 
dashed line refers to the present BR32 or BR31 experimental limit, respectively l29ll . It is clear from these two plots 
that the ongoing MEG experiment stands the best chance of confirming the predictions for or eliminating the GK and 
AB models. Even with a super-B factory, the present experimental bounds on the BR32 and BR31 branching ratios 
can only be lowered by one or two orders of magnitude at most 13011 . 

The above figures apply in the case where Aq = is selected for the common trilinear scalar coupling. We show in 
rt25ll that as |Ao| is allowed to depart from zero, the predicted branching ratios increase. Hence the line segments in 
Figs. 4 and 5 represent lower limits and extend upward somewhat at 45° as \Aq\ is increased. 



8. LEPTON FLAVOR VIOLATION IN jU - e CONVERSION 

Lepton flavor violation can also occur in the fl — e conversion process in 81 Ti, where fl + Ti—* e + Ti. The \i — e 
conversion branching ratio is the conversion rate scaled by the capture rate for the process, /1 + 77 — ► + Sc. The 
one-loop diagrams involving 7, Z and Higgs penguins all contribute along with box diagrams, but in the CMSSM 
scenario the / penguin has been shown to dominate [31]. We show two such diagrams involving slepton-neutralino 
and sneutrino-chargino loops in Fig. 6, where the effects of the virtual N[ and N[ with their Yukawa couplings appear 
in slepton loops. 

The n — e conversion branching ratio on Ti vs. BR21 is plotted in Fig. 7 for the five GUT models, where the tighter 
WMAP dark matter constraints have been imposed, again for the case of Aq = 0. The present conversion branching 
ratio limit for the Ti experiment IB2I1 is shown at 4 x 1CU 12 . It is projected that such an experiment will be able to 
reach down to 10~ 17 in a first round and down to 10 -18 in a second generation experiment. While the expectation is 
that the MEG experiment may be able to reach a limit of 10~ 13 branching ratio for )X — > ey, it is apparent that a fi — e 
conversion experiment which eventually reaches a limit of 10~ 18 will be considerably more powerful. Potentially such 
an experiment would be able to eliminate all five models, and probably all such SUSY GUT models in the CMSSM 
scenario, if no positive signal is found. 




FIGURE 6. Examples of Feynman diagrams for slepton - neutralino and sneutrino - chargino contributions to fl — e conversion 
in SUSY models with slepton mass insertions. 
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FIGURE 7. Branching ratio predictions for fl — e conversion vs. branching ratio predictions for fl — » e + y in the five models 
considered. The more restrictive WMAP dark matter constraints apply for the thick line segments shown. Note that the predictions 
for the CM and CY models nearly overlap. 



9. CONCLUSIONS 

We have tried to differentiate models based on sin 2 and charged lepton flavor violation predictions. Our study was 
initially based on 63 models available in the literature prior to June 2006. There we found that a normal neutrino 
mass hierarchy is preferred 3:1. Moreover, future Double CHOOZ and Day a Bay reactor experiments will be able to 
eliminate roughly half of the 63 neutrino models surveyed, if their sensitivity reaches sin 2 2 613 ~ 0.01 as planned. Still 
of the order of five models have similar values for sin 2 613 in the interval 0.001 - 0.08. We have suggested that charged 
lepton flavor violation experiments may be able to further distinguish them. If the now-running MEG experiment sees 
positive signals for pi — > ey, all non-SUSY models or models which do not involve new physics will be ruled out. 

We then narrowed our study to five predictive 5(9(10) SUSY GUT models in the literature. All five models have 
type I seesaw mechanisms implying normal hierarchy. Their predictions for sin 2 20b lie in the ranges of 0.05 for the 
Chen-Mahanthappa model, 0.01 for the Albright-Barr, Cai-Yu, and Dermisek-Raby models, and 0.001 for the Grimus- 
Kubock model. Previous studies of generic SO(10) models have concluded that the LFV branching ratios depend 
critically on ©13 and the heaviest right-handed neutrino mass, M3. Here we find that M3 appears to be more important. 
If the MEG experiment can reach an upper bound of 10~ 13 for the — > ey branching ratio, it will rule out the Grimus- 
Kubock and Albright-Barr models if no positive signal is seen. If a jx — e conversion experiment can be performed and 
reach a branching ratio limit of 10~ 18 as projected, it can potentially rule out all five models considered. 



ACKNOWLEDGMENTS 

The work reported on here was carried out in collaboration with Mu-Chun Chen. The author thanks the members of 
the Theory Group at Fermilab for their kind hospitality. Fermilab is operated by the Fermi Research Alliance under 
contract No. DE-AC02-07CH1 1359 with the U.S. Department of Energy. 



REFERENCES 



1. M. Maltoni, T. Schwetz, M.A. Tortola and J.W.F. Valle, New J. Phys. 6, 122 (2004), |arXiv:hep-ph/0405172 updated 28 Sep 
2007]. 

2. P.F. Harrison, D.H. Perkins, and W.G. Scott, Phys. Lett. B 530, 167 (2002). 

3. D.N. Spergel et al. [WMAP Collab.], Astrophys. J. Suppl. 170, 377 (2007). 

4. P. Minkowski, Phys. Lett. B 67, 421 (1977); T. Yanagida, in Proceedings of the Workshop on the Unified Theory and Baryon 
Number in the Universe, Tsukuba, Japan 1979, edited by O. Sawada and A. Sugamoto (KEK, Tsukuba, 1979), p. 95; M. 
Gell-Mann, P. Ramond, and R. Slansky, in Supergravity, (P. van Nieuwenhuizen and D. Z. Fredman, eds.), North-Holland, 
Amsterdam, 1979, p. 315; S.L. Glashow, in Proceedings of the 1979 Cargese Summer Institute on Quarks and Leptons, (M. 
Levy, J.-L. Basdevant, D. Speiser, J. Weyers, R. Gastmans, and M. Jacob, eds.), Plenum Press, New York, 1980, p. 687; R.N. 
Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

5. R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980); J. Schechter and J.W.F. Valle, Phys. Rev. D 22, 2227 
(1980). 

6. B. Pontecorvo, Sov. Phys. JETP 6, 429 (1957) [Zh. Eksp. Teor. Fiz. 33, 549 (1957)]; Z. Maki, M. Nakagawa and S. Sakata, 
Prog. Theor. Phys. 28, 870 (1962). 

7. For recent reviews, see, e.g. M.-C. Chen and K. T. Mahanthappa, Int. J. Mod. Phys. A 18, 5819 (2003); AIP Conf. Proc. 721, 
269 (2004); G. Altarelli and F. Feruglio, New J. Phys. 6, 106 (2004); S.F. King, Rept. Prog. Phys. 67, 107 (2004); Z. z. Xing, 
Int. J. Mod. Phys. A 19, 1 (2004); R.N.Mohapatra, New J. Phys. 6, 82 (2004); A. Yu. Smirnov, Int. J. Mod. Phys. A 19 1180 
(2004); R.N. Mohapatra and A .Y. Smirnov, Ann. Rev. of Nucl. and Part. Sci., 56, 569 (2006); R.N. Mohapatra et al, "Theory 
of Neutrinos: A White Paper," |hep-ph/05 1 0213] ; A. Strumia and F. Vissani, |hep-ph /0606054 1. 

8. C.H. Albright and M.-C. Chen, Phys. Rev. D 74, 1 13006 (2006). 

9. SA. Dazley [Double CHOOZ Collab.], in Proceedings of NuFACT05, Nucl. Phys. Suppl. 155, 231 (2006); J. Cao [Daya Bay 
Collab.], ibid. 155, 229 (2006). 

10. A. Baldini, Nucl. Phys. B Suppl. 168, 334 (2007). 

11. M.L. Brooks et al. [MEGA Collaboration], Phys. Rev. Lett. 83, 1521 (1999). 

12. B.W. Lee and R.E. Shrock, Phys. Rev. D 16, 1444 (1977). 

13. For references to the CMSSM, cf. K.A. Olive, submitted to the SUSY07 proceedings | arXiv: 0709. 33031 . 

14. F. Borzumati and A. Masiero, Phys. Rev. Lett. 57, 961 (1986). 

15. S.T. Petcov, S. Profumo, Y. Takanishi, and C.E. Yaguna, Nucl. Phys. B 676, 453 (2004). 

16. J. Hisano, T. Moroi, K. Tobe, and M. Yamaguchi, Phys. Rev. 53, 2442 (1996). 

17. A. Masiero, S.K. Vempati, and O. Vives, Nucl. Phys. B 649, 189 (2003); New J. Phys. 6, 202 (2004); F. Deppisch, H. Pas, A. 
Redelbach, R. Ruckl, and Y. Shimizu, Eur. Phys. J. C 28, 365 (2003); F. Deppisch, H. Pas, R. Ruckl, and A. Redelbach, Phys. 
Rev.D 73, 033004 (2006); E. Arganda and M.J. Herrero, Phys. Rev. D 73, 055003 (2006); L. Calibbi, A. Faccia, A. Masiero, 
and S.K. Vempati, JHEP 7, 012 (2007). 

18. J.A. Casas and A. Ibarra, Nucl. Phys. B 618, 171 (2001). 

19. C.H. Albright and S.M. Barr, Phys.Rev. D 64, 073010 (2001); C.H. Albright, Phys. Rev. D 72, 013001 (2005); Erratum-ibid.D 
74, 039903 (2006). 

20. M.-C. Chen and K.T. Mahanthappa, Phys. Rev. D 70, 1 13013 (2004). 

21. Y. Cai and H.-B. Yu, Phys. Rev. D 74, 1 15005 (2006). 

22. R. Dermisek and S. Raby, Phys.Lett. B 622, 327 (2005). 

23. W. Grimus and H. Kiihbock, Phys. Lett. B 643, 182 (2006). 

24. A. Pilaftsis and T.E.J. Underwoo d, Nucl. Phys. B 692, 303 (2004). 

25. C.H. Albright and M.-C. Chen, larXiv:0802.42281 . 

26. Review of Particle Physics, W.-M. Yao et al., Jour, of Phys. G 33 1 (2006). 

27. P. Binetruy, G. Girard'i, and P. Salati, Nucl. Phys. B 237, 285 (1984); K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991). 

28. L.S. Stark, P. Hafliger, A. Biland, and F. Pauss, JHEP 0508, 059 (2005). 

29. B. Aubert et al. [BABAR Collab.], Phys. Rev. Lett. 96, 041801 (2006); ibid. 95, 041802 (2005). 

30. A.G Akeroyd et al. [SuperKEKB Physics Working Group], |hep-ex/0406071 ]. 

31. E. Arganda, M.J. Herrero, and A.M. Teixeira, larXiv:0707.29551 . 

32. C. Dohmen et al. [SINDRUM II], Phys. Lett. B 317, 631 (1993). 



